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Abstract 

A set of one dimensional interfaces involving attachment and detachment of 
fc-particle neighbors is studied numerically using both large scale simulations and 
finite size scaling analysis. A labeling algorithm introduced by Barma and Dhar 
in related spin Hamiltonians enables to characterize the asymptotic behavior of the 
interface width according to the initial state of the substrate. For equal deposition- 
evaporation probability rates it is found that in most cases the initial conditions 
induce regimes of saturated width. In turn, scaling exponents obtained for initially 
flat interfaces indicate power law growths which depend on k. In contrast, for 
unequal probability rates the interface width exhibits a logarithmic growth for all 
k > 1 regardless of the initial state of the substrate. 



PACS numbers: 05.40.-a, 05.70.Ln, 64.60.Ht 

KEY WORDS: irreducible string; saturated regimes; logarithmic growth. 
Published in J. Stat. Phys. 103, 395 (2001). 



1 



Studies of interface growth have acquired a major momentum in the past decade 
[|]. 0, |3|] , finding a rich variety of applications. These range from molecular beam epi- 
taxy j|] to magnetic flux lines in type-II superconductors [1], || and polymer physics 
||. Most numerical analysis and theoretical investigations on growing structures 
have emphasized the onset of scaling regimes which emerge at both large time and 
length scales. This enables a classification of nonequilibrium interfaces by universal- 
ity classes characterized by a set of scaling exponents which, independently of the 
microscopic growth rules, dominate the late evolution stages of these processes. 

However, the issue of universality classes in these problems has been steadily 
brought into question, particularly in higher dimensions |7|, ||, ||, pL0| . Since no 
systematic analytical method is yet available to understand large scale fluctuations 
in nonequilibrium systems, two main research lines have been addressed to probe 
the universality hypothesis. They involve either the simulation of discrete models 
[PH , |ll|, [12], [13], [14j] or the numerical solution of phenomenological Langevin-type 
equations of growth driven by Gaussian noise |1| . This is the case of the ubiquitous 
Kardar-Parisi- Zhang (KPZ) equation [17| which has been applied to characterize a 
vast body of far from equilibrium processes |1|, [2], [!| . 

In this letter we rather follow the first category of investigations and introduce a 
set of restricted solid-on-solid (RSOS) growth models in one dimension entailing an 
exceptionally large number of conservation laws [|18|]. By mapping the problem onto 
a quantum spin system, it will turn out that for equilibrium situations, not only 
each model defines a "class" of its own but furthermore, their asymptotic properties 
become extremely sensitive to slight deviations from initially flat substrates. In 
contrast, non-equilibrium situations are rather robust but exhibit a much slower 
dynamics. This will mark a clear difference between fluctuations in equilibrium and 
non-equilibrium cases. 

The basic process considered is a simple yet nontrivial extension of ballistic ag- 



gregation models |14j| . Here, we include both deposition and evaporation of extended 
objects such as dimers, trimers, ... fc-mers, at random locations of a one-dimensional 



interface As usual, the state of the system at a given time is defined by a set 
of single valued functions {hi(t), ... , h^t)}, where h n is the height of the interface 
measured from a reference level at site n & {1,2, ... , L}. To prevent the divergence of 
surface fluctuations at large times we also impose a RSOS constraint throughout the 
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process namely, | h n+ i — h n \ = 1 , Vt . Therefore, a fc-mer deposition (evaporation) 
attempt with rate e (e ; ) is successful provided a plateau with at least ^-successive 
local minima (maxima) is selected from the interface. Fig. 1 illustrates these mi- 
croscopic rules for the case of dimers. Notice that evaporation apply whether or 
not the targeted fc-adjacent maxima were created together, a feature which allows 
for reconstitution of fc-mers. Also, it is helpful pointing out that the dynamics of 
these interfaces is equivalent to that defined in the 2fc-mer deposition-evaporation 
problem studied in Ref . , by flipping the spins on all even sites, say, of this latter 
system. 

In studying the scaling regimes of growing structures it is useful to consider the 
mean square fluctuations W 2 of the average interface height (h(t)) , i.e. W 2 (L,t) = 
t J2j( [hj{t) — ( h(t) ) ] 2 ) , which yields a measure of the interface width. From the 
assumption of a scaling form of height-height correlation functions along with a 
single time-dependent length scale ~ t 1 ^ for growth starting from an initially flat 



interface, it can be argued |L9| that W scales as 

W(L,t) = L<f(t/L z ), (1) 
where the scaling function /(c) satisfies 



/(c) „jc</« for c<<1 (2) 
I const for c^> 1 . 

It follows immediately that finite systems saturate as W oc , whereas in the 
thermodynamic limit the asymptotic growth is dominated by the exponent (3 = (/z, 
that is W oc t 13 . The exponent ( describes the roughness dependence of the surface 
width on the typical substrate size. In turn the exponent z, often known as the 
dynamic exponent, gives the fundamental scaling between length and time. 

However, already for k > 2 it will turn out that there are initial substrates whose 
asymptotic widths neither follow a power law nor exhibit the universal characteristics 
embodied in Eqs. (0) and (0). As we shall see, the investigation of these issues 
becomes rather systematic in the language of the spin systems discussed below. 

a. Spin representation — We now exploit the well known mapping between 
RSOS interface dynamics and quantum spin systems [Q. Associating the height 
difference s n = h n+ i — h n to an eigenvalue of the ^-component, say, of the Pauli 
operator a n for site n, then all intrinsic properties of the interface can be analyzed in 
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terms of direct products \s) = ... , sl) of 1/2 spinors. By construction it is clear 
that up to a constant, the interface heights are obtained as h n = hi+Y^=\ Sj , where 
hi — hi = s\ upon imposing periodic boundary conditions (PBC). Note that the 
latter force a vanishing total magnetization throughout the underlying spin kinetics. 

Introducing spin-1/2 raising and lowering operators <r+, a~ , the stochastic evo- 
lution of our fc-mer interfaces at subsequent times can be calculated from the action 
e -Ht Q f quantum "Hamiltonian" 

H = EM + ^H^+A- 1 ) > ( 3 ) 

3 
k 

4 = IK+2i-2^ + 2i-l- (4) 



i=l 



Here, adsorption (desorption) of fc-mer plateaus with probability e (e') is described 
by the effect of the A> (A ) operators. Conservation of probability requires in turn the 
appearance of 2/c-field correlators, already diagonal in the a z or particle (monomer) 
representation. We address the reader to Ref. [f20| for a more detailed derivation in 
related systems. 

b. Conservation laws — The advantage of the spin representation is that one 
can identify a large number of subspaces which are mutually disconnected by the 



if-dynamics. Following Ref. ||18|| , we define accordingly a reduction rule by looking 
for the occurrence of groups of 2fc-opposite consecutive spins in a given state \s). 
Each occurrence, if any, is deleted so the length of the remaining spin configuration 
is reduced in 2k bits per deletion. This procedure is applied recursively, until we are 
left with a string that cannot be further reduced. Such an object corresponds to an 
irreducible string I\s) jTS]. For example, in the dimer case we have I\ TTIItllT) = 
TTilTliT (irreducible block of Fig. 1), i\ TITITUT) = TUT , /| TITITITI) = 0, 
the null string. 

The key observation is that so long as e, e' > 0, two spin configurations \s) 



and \s') belong to the same subspace <^=^> I\s) = I\s') [0]. Thus, treated as 
a dynamical variable, the irreducible string is a constant of motion and provides 
a unique label for each i7-invariant subspace. To illustrate this point consider a 
sector labeled by xi > ■•• ; Xc wr th £ < L irreducible characters x — ±1 • Let {xj(t)} 
denote their positions in the full L-site lattice before applying the reduction rule. 
The interface dynamics then induces a random walk of these characters such that 
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Xj +1 (t) > Xj(t) Vt . Since H changes Xj by multiples of 2k lattice spacings, the spin 
array between sites Xj and Xj+i can be ultimately reduced to a null string. Thus, 
we see that {xj(t)} constitute the slow modes of our fluctuating interfaces whereas 
their characters remain unaltered throughout the dynamics. The above reduction 
algorithm also enables to count the total number of invariant sectors, both jammed 
(strings which are already irreducible), and unjammed. A straightforward analysis 
given as in [ffH] shows that these sectors increase as fast as X L for large L, where A 



is the largest root of A = 2 A — 1 . 

c. Saturated regimes — It is then natural to ask whether the dynamics is affected 
by this rather unusual partition of the full phase space. At least in the simpler case, 



when e = e' |[21|| , it can be argued that W remains bounded whenever we start from 
a substrate of finite width and the length £ of the corresponding irreducible string 
is a finite fraction of the lattice size. Specifically, consider two arbitrary points A 
and B in a ring of L sites. The most general configuration of spins determining the 
height difference h& — Hb can be denoted as 



X 1 Wx x 2 w 2 ... X n _i W n _i x n 



U2 



where Xi are the positions of the irreducible characters between A and B, wi are 
spin patches reducible to null strings, and U\ (u 2 ) is the part of the irreducible string 
between the previous (next) irreducible character x (x n+ i ) that lies to the left of 
A (right of B). To calculate the height difference between A and B, we can ignore 
all Wi patches as they necessarily have equal number of up and down steps. On the 
other hand, if the initial substrate has bounded heights, the sequence Xi(t) ... x n (t) 
can only give a bounded contribution to the height difference Vt. The remaining 
contributions at subsequent times comes from u\ and u 2 . If u\ is part of a reducible 
substring of length 2m (between xq and x\), the maximum height difference it can 
have is at most m. Also, for e — e' the detailed balance condition holds in the steady 
state. For such equilibrium situation it has been shown that m has an exponential 
distribution, provided that the density C/L of irreducible characters is held finite 
[ |l"8|j . Consequently, u\ yields only a finite variance. Similar considerations apply to 
u 2 . Hence, the height difference between two arbitrary points has bounded variance 
Vt and L. Therefore, we conclude that the width of the substrate remains finite, 
or equivalently, in the steady state W should exhibit a saturated regime. However, 
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notice that if the irreducible string is very short, i.e. C/L — > , fluctuations in the 
(xj+i —Xj) distances become large and the above reasoning for u±, u<i does not hold. 

The case k = 2 is typical and already exhibits such sector dependent behavior. 
Thus, we consider only this case in detail. 

d. Numerical results — Armed with this labeling methodology along with the 
spin-height mapping discussed above, we studied the dimer interface width W in a 
number of representative subspaces. Fig. 2 displays the width behavior obtained 
from Monte Carlo simulations using five initial scenarios. These were embedded 
homogeneously on the Neel state of 2.1 x 10 5 lattice sites with PBC, and corre- 
spond to (1) the null string, i.e. the usual flat substrate; (2) [ttll] i ^ 60 that is, the 
irreducible string constructed by repeating the unit cell L/60 times; (3) [Tt!!] L ^ 8 ; 
(4) [Ttlllt] L ^ 14 ; (5) a random configuration formed by deposition-evaporation of 
monomers on the Neel state. Setting e = e', the averages were taken over 100 
histories of 10 5 steps per height. 

Although for e = e' the interface of sector (1) moves on the average with zero 
velocity [see lower inset of Fig. 3(a)], alike the monomer case (k = 1) it can 
grow to arbitrarily large heights. In fact, this is consistent with our data which 
support a power law growth W 2 oc t 2f3 extended over more than two decades with 
an exponent (3 ~ 0.12(1) . On the other hand, sectors (3) and (4) clearly exhibit 
the saturated regimes conjectured above. They fit well with an exponential form 
W 2 oc exp (— t~ a ) , as is shown in the semi-log inset of Fig. 2 upon setting a ~ 1/4 
and 1/3 respectively in sharp contrast with the expectations of Eqs.(§]) and @. 
Thus, we see that certain deviations from a perfectly flat substrate turn into drastic 
changes at large times. As it was referred to above, the latter depend on both the 
density p = C/L and particular sequence of irreducible characters considered p3fl . 

The situation of sectors (2) and (5) is less clear. Here p is still finite (though 
much smaller than the irreducible density of the previous non null sectors), and 
consequently a saturated regime would be expectable. However, also here the data 
extend nearly over two decades following power law growths with (3 ~ 0.07 and 
0.04 respectively (dashed slopes of Fig. 2). One may attempt to interpret these 
results in terms of a temporal crossover, although to reach an asymptotic saturated 
regime such a crossover should be exceptionally slow. A numerical test in favor of 
this latter possibility is provided by the average interface height ( h(t) ). Notice that 
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given a distribution of heights P(h\, ... , h L , t), the mean height can be evaluated as 

dimer dimer 



dt(h) = -E eP (- fn, h£) hi , hf> ... , t) - e'P(... hf] hf] hf] hf] hP ... , t) , (5) 

L \h),i 

where hf — hi+j , and the overbraces involve the heights which form a dimer, say 
hi_i, ... ,/ij +3 (see Fig. 1). For e = e', it follows that (h(t)) must saturate to a 
finite value at large times, as all accessible height configurations are equally likely 



in the steady state, irrespective of the sector considered [ 21] . In fact, Fig. 3 (a 



displays saturation regimes in most sectors. However, it also reveals the existence of 
a rather huge temporal crossover in sector (2) (more than 10 5 steps), along with a 
medium one in sector (5) (about 10 3 steps), possibly responsible for the unexpected 
width behavior observed in Fig. 2. It is important to note that for e > e' the 
average interface height behaves quite differently, at least within our accessible time 
scales. The latter can grow indefinitely with a non zero velocity, as suggested by 
the asymptotic logarithmic increase of (h(t)) observed in Fig. 3 (b) for all sectors. 

Before discussing such logarithmic behavior (see section e below), we complete 
our study of the case e = e' and analyze the late growing stages of higher k-mer 
interfaces, now focusing attention on initially null strings, the common flat substrate. 
After averaging over 100 histories our results indicate a slight dependence of j3 on 
the value of k (however, see below the variation of z), which nonetheless can be 
clearly distinguished over five step decades, as indicated by the slopes of Fig. 4. 
The measured values of (3 are presented in Table I for k = 2, 4 and 8. 

To enable a complementary check of these exponents we turn to a finite size scal- 
ing analysis of the interface width. The reader's attention is directed to Figs. 5(a)-(c) 
where we display in turn the results of the simulations carried out for k = 2, 4 and 
8. Here, the averages were taken over 2 x 10 4 samples employing periodic chains 
of L = 400, 640 and 800 sites. Using the estimates of the dynamic and rough- 
ening exponents {z, (} given in Table I, we obtained a fair data collapse towards 
the (squared) phenomenological scaling functions conjectured by Eqs. (|1|) and (|2|) 
[ |I~9|j . As expected, the resulting values of (/z are in reasonable agreement with the 
previous /3-set obtained in the thermodynamic limit. So we finally see that as long 
as the deposition-evaporation rates are held equal, the emerging scaling exponents 
are in fact nonstandard and give rise to rather unusual dynamic length scales t l l z 
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which decrease markedly with k. 

e. Logarithmic growth, e / t' — The logarithmic dependence of the average 
interface height (h(t)) obtained in Fig. 3(b) at large times, may be understood 
by means of the following heuristic considerations. For clarity of argument, assume 
e 3> e'. The random deposition leaves some single spaces between segments occupied 
by dimers. This single spaces cannot be filled by subsequent deposition of dimers, 
with the result that the height profile after many depositions resembles a "moun- 
tain landscape" where no further depositions are possible. Such profile is clearly 
evidenced by the evolution snapshot given in the upper panel of Fig. 6. 

If e' — 0, this configuration cannot evolve any further, and the configuration is 
stuck. If e' is non-zero, one can evaporate away from one side of the mountain, all 
the way to the base, and then redeposit the base by dimer shifting it one unit (right 
or left, whatever is allowed). Then, rebuilding on this side, one gets a different 
profile (dotted line in the lower panel of Fig. 6). This rebuilding process leads to 
sometimes extinction of small mountains (which merge with bigger ones), and thus 
the average mountain size M(t) increases with time. The time dependence of M 
can be estimated as follows. For a single rearrangement event described above, one 
needs order M evaporation events. Hence by standard activation arguments, this 
rate is (e'/e) M . One needs order M 2 such events for total width to change by M 
(since change occurs either way). Hence, the time for the average mountain size to 
change from M to 2M is of order M 2 (e/e') M , so we get t(2M) ~ M 2 {e/e') M . On 
the other hand, from simple geometry both the average mountain height and the 
width W scale as M (i.e. the roughness exponent C is 1). Therefore, this yields 

W{t) ~ log(t) , e^e'. (6) 

These conclusions can be extended for all substrates (whatever the sector or irre- 
ducible string considered), as well as V/c > 1. In fact, the semi- log data displayed in 
Fig. 7 corroborates the above conjecture for k = 2, 4, and 8 starting from empty sub- 
strates. Also, the upper inset confirms this logarithmic growth for k = 2 in sectors 
(2)-(5). A similar one-dimensional model showing logarithmically slow coarsening 
was studied in . 



To summarize, we have studied a family of simple growth models yet exhibiting 
highly nontrivial behavior. For e = e' our finite size analysis indicates that each 
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k-mer process might be regarded as a growth class of its own dominated by different 
dynamic exponents z(k). The large systems explored within the time scales of 
Fig. 4 yielded growth exponents consistent with the (3 values obtained from finite 
substrates (Fig. 5). The concept of irreducible string ]TB| played a crucial role in 
determining the saturated behavior entailed by the initial conditions (Fig. 2), as 
opposed to the much simpler dynamics of monomer interfaces. However, the late 
dynamic stages of small strings (C <C L) involve slow temporal crossovers for which 
further numerical efforts will be required to clarify such situations. Generalizations 
of this labeling algorithm to higher dimensions are clearly desirable though its mere 
existence appears difficult to elucidate. Whether or not an exponential number of 
conservation laws show up in d > 1 yet remains an open question. 

For the non-equilibrium situation (e ^ e'), Fig. 7. indicates a rather robust 
dynamics. Here, the interface width grows logarithmically VA; > 1 irrespective 
of the initial state of the substrate. In view of these drastic differences between 
fluctuations in equilibrium and non-equilibrium cases, the transient dynamics of 
regimes with e ~ e' still needs further investigations. 

It is a pleasure to thank R.B. Stinchcombe, M. Barma and T.J. Newman for help- 
ful comments. Also, the remarkable observations made by the referee contributed 
to improve this work. The author acknowledges support of CONICET, Argentina. 
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k 


P (L -> oo) 


z 


c 


2 


0.12(1) 


2.5(10) 


0.3(1) 


4 


0.10(1) 


2.8(10) 


0.3(1) 


8 


0.09(1) 


3.3(10) 


0.3(1) 



Table 1: Growth exponent (3 measured in large /c-mer substrates (Fig. 4), along 
with finite size scaling estimations of dynamic and roughening exponents z and ( 
(Fig. 5), obtained for equal deposition-evaporation rates. 




Figure 1: Schematic representation of dimer deposition- evaporation onto a RSOS 
interface. The equivalent spin-| (s n = h n+ i — h n ), or hard core particle dynamics 
involves two left (right) particle hopping at a time with rate e (e') for dimer adsorp- 
tion (desorption). The array between brackets illustrates an irreducible string (see 
text). 
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Figure 2: Evolution of W 2 in subspaces (l)-(5) defined in text. Here, they correspond 
to (1) circles (null string); (2) romboids; (3) triangles; (4) pentagons and; (5) squares. 
The slopes of upper straight lines indicate a power-law asymptotic behavior t 213 with 
(3 ~ 0.12 (top), 0.07 (middle), 0.04 (bottom). These two latter cases however, might 
be undergoing a slow temporal crossover (see text). In turn, subspaces (3) and 
(4) fit appropriately an exponential form W 2 oc exp (— t~ a ) with a ~ 1/4 and 1/3 
respectively, as shown in the inset. 
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Figure 3: Average interface height (h) for k = 2 with (a) e — e', measured from the 
flat interface level [(h(0)) = 0.5]. Sectors (1) (bottom curve of lower inset), (4) and 
(3) (from top to bottom in the upper inset), display a fast relaxation toward the 
expected saturated regime. In contrast, sectors (2) (romboids) and (5) (top curve 
of lower inset), undergo a slow temporal crossover, (b) Average height {h) for k = 2 
with e'/e = 0.4 . All sectors exhibit logarithmic growth at large times as opposed to 
the saturated regimes in (a). Solid lines are guide to the eye whereas all symbols 
are taken as in Fig. 2. 
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2 3 4 5 
1 10 10 10 10 10 
Step 



Figure 4: Growth of interface width W 2 using different A;-mers, L = 2.1 x 10 5 lattice 
sites and e — e'. For comparison, the uppermost curve shows the monomer case 
whereas lower curves correspond to k — 2 , 4, and 8, in descending order. For k > 1, 
straight lines are fitted with the values of 2(3 given in Tab. I 
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Figure 5: Finite size scaling of interface width W(L,t) using null strings and e — e' 
for (a) k — 2, (b) k = 4 and, (c) k = 8. Sizes L = 400, 640, and, 800 are denoted 
respectively by squares, circles and triangles. The data collapse was attained upon 
tuning the set of exponents {z, (} given in Table I. The straight lines are fitted with 
the values of 2f3 obtained in the thermodynamic limit. 
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Figure 6: Evolution snapshot for dimers after 10 4 steps for e'/e = 0.4. Only 500 sites 
are shown (upper panel). The dotted lines in the schematic representation denotes 
the rebuilding process described in the text (lower panel). 
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Figure 7: Logarithmic growth of the interface width for e'/e = 0.4 starting from a 
flat substrate. Straight lines over circles (k — 2), romboids (k — 4), and triangles 
(k = 8) are guides to the eye. The upper inset displays similar results for k = 2 in 
sectors (3) (triangles), and (5) (squares), which are closely and respectively followed 
by those of sectors (4) and (2) (solid lines). 
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